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ABSTRACT 

A significant fraction of the stellar population in the cusp around central black holes of 
galaxies consists of compact remnants of evolved stars, such as white dwarfs, neutron 
stars and stellar mass black holes. We estimate the rate of capture of compact objects 
by massive central black holes, assuming most spiral galaxies have a central black 
hole of modest mass (~ IO^Mq), and a cuspy spheroid. It is likely that the total 
capture rate is dominated by nucleated spirals. We estimate the flux of gravitational 
wave radiation from such coalescences, and the estimated detectable source count for 
proposed space-based gravitational wave observatories such as LISA. About one event 
per year should be detectable within 1 Gpc, given very conservative estimates of the 
black hole masses and central galactic densities. We expect 10 2 -10 3 detectable sources 
at lower frequencies (10 -4 Hz) "en route" to capture. If stellar mass black holes are 
ubiquitous, the signal may be dominated by stellar mass black holes coalescing with 
massive black holes. The rate of white dwarf-white dwarf mergers in the cores of 
nucleated spirals is estimated at ~ 10 -6 per year per galaxy. 

Key words: galaxies-nuclei, black holes: stellar-mergers 



1 INTRODUCTION 

There is increasingly strong evidence that massive black 
holes are found in the centers of galaxies (see eg. Kormendy 
& Richstone 1995, Tremaine 1995, Rees 1990). Observa- 
tional evidence, and theoretical considerations, indicate that 
the masses, M/j, range from under 10 6 Mq to over 10 9 Mq. 
Black holes are thought to form during, or possibly before 
the formation of the host galaxy. Early infall of low angu- 
lar momentum material, and angular momentum transport 
in the disk of gas surrounding the proto-black hole fuels 
the early stages of black hole growth, with accretion rates 
^ lM y _1 inferred (see eg. Rees 1990). The lifet ime of 
the resultant quasar is uncertain, it is possible that less than 
1% of galaxies host quasars, with lifetimes ^ 10 9 y, but more 
likely that most galaxies undergo shorter periods of activity 
(see eg. Haehnelt and Rees 1993). 

If the QSO duty cycle is low, then most galaxies must 
have undergone moderately short periods of accretion onto 
a central black hole, in order to account for the total QSO 
numbers, and a typical galaxy has a ^ 10 7 Mq cen- 
tral black hole. If the QSO/AGN duty cycle is high, then 
remnants in active galaxies are mostly super-massive and 
as many as ^ 99% of (non-active) galaxies will have moder- 
ate mass (M/j ~ 1O 6 M0) black holes that did not undergo 
prolonged strong accretion episodes. It is interesting to note 
that for our own galaxy and the nearby M32 dwarf elliptical, 
the best estimates of the mass of any central black hole is 
^ 3 x 1O 6 M (Bender et al 1996, van der Marel et al 1996, 



Eckart and Genzel 1996). The true black hole mass function 
probably lies between the extremes mentioned, with a tail 
of very high mass black holes in a few galaxies, and some 
unknown distribution of masses in normal galaxies. A useful 
guess for the black hole mass function assumes the mass of a 
central black hole is some fraction of the mass of the stellar 
spheroid, = fgL, where f g PS 1O -2 ' 5 M0/L0 , and L is 
the galaxy luminosity (Tremaine 1995). We assume for now 
that M/j cx L and a Schecter (1976) luminosity function 

$(L) = $o(^)V(^)f, (1) 

where U = 1.8xlO lo L , $ = 0.008 Mpc" 3 and a = -1.07 
(Tremaine 1995, Efstathiou et al. 1988). 

In general a dense cusp of stars (and possibly gas) forms 
around any central mass in a galaxy (Peebles 1972, Young 
1980). Stars in these cusps will occasionally come close 
enough to the black hole to collide with each other or be 
swallowed by the black hole (see eg Frank and Rees 1976, 
Hills 1975). For black holes with masses near the low end of 
the range expected (M/j ~ 1O 6-7 M0) stars and (sub)giants 
are tidally disrupted if they come close to the black hole; 
more massive black holes may swallow stars whole. We ex- 
pect ~ 10% of stars in an evolved population, such as is 
observed in elliptical galaxies and the bulges of spiral galax- 
ies to be white dwarfs(WDs), another <J 1% of stars may 
be neutron stars(NSs) or even low mass (7 — 100 Mq) black 
holes(LMBHs), depending on the initial mass function and 
fate of evolved massive stars, which is uncertain. These 
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evolved stars will also encounter the central black hole, but 
won't be disrupted, and may be promising sources for grav- 
itational waves detectable by proposed detectors such as 
LISA (Danzmann et al. 1993, Hough et al. 1995). 

1.1 Stellar encounters with central black 
holes 

A star, radius r* , mass to*, would be tidally disrupted 
by a black hole mass within a radius rj< PS 2 x 

(M ft /m*) 1/3 r*. (Rees 1988, Evans and Kochanek 1989, 
Frolov et al. 1994). For a main sequence star of solar 
mass, encountering a 10 6 Mg Mq black hole, PS 1.4 x 
cm(fi 50rsM e 2/,S . In general the central region 
of a galaxy will contain both main sequence and evolved 
stars, with a range of masses and radii. We need to assume 
some global initial mass function for the stars, here taken 
to be the Salpeter mass function, dN*/dm* = m~^ 1 ~ x * , 
where x* = 1.35. A typical galaxy will contain an evolved 
population, with a turnoff mass for the main sequence of 
0.7 — I.OMq, plus possibly an additional younger popula- 
tion. The evolved population will typically consist of about 
0.2% neutron stars of 1.4Mq, about 7% white dwarfs, with 
masses from 0.5-1.3Mq, and possibly about 0.03% few Mq 
black holes. As we discuss later, however, dynamical effects 
may alter these proportions within the central cusp. 

A degenerate, dark, compact remnant of the same mass 
as a main sequence star has a much smaller radius. For white 
dwarfs, r\yjj ss 0.01 Rq, for neutron stars, rjys ~ 10 6 cm, 
much less than r* . Consequently the white dwarf tidal ra- 
dius is smaller than that for a main sequence star of the same 
mass by approximately ryyD l r * ■ Since the Schwarzschild 
radius of a black hole, r s = 3 x lO n M ft /(lO 6 M )cm, in- 
creases linearly with , we see that for a black hole with 
Mh 1O 6 M0, rs > tt f° r white dwarfs. For the values of 
M/j considered here, neutron stars are not subject to tidal 
disruption; nor, obviously, are any low mass black holes, 
which may also encounter the central black hole. Compact 
stellar remnants can orbit around galactic black holes inside 
the tidal disruption radius for stars. 

The orbital period for a solar type star around black 
hole at is = 10 4 s, independent of M^. The charac- 
teristic frequency, f c , of gravitational radiation associated 
with the orbit of compact remnants about the black hole 
inside is then 10 -4 -10 -2 Hz. Gravitational radiation of 
this frequency is not readily detectable by Earth bound de- 
tectors but is well matched to the expected sensitivity of 
space based detectors (Thorne 1987, 1995, Danzmann et al. 
1993, Hough et al. 1995, Haehnelt 1994). The sensitivity of 
these detectors cuts off below ~ 10 -3 Hz, and the gravita- 
tional radiation rate by a star in the innermost stable orbit 
around a massive black hole scales as M^' 2 " . For this reason 
we are primarily interested in moderate mass central black 
holes, 0.5 <, M 6 <, 5. 

The corresponding characteristic amplitude of the grav- 
itational waves, h c , for a source at distance dQ pc = d/l Gpc, 
M h > to*, is 

where P is the orbital period in seconds and g is a geometric 
factor of order unity (Thorne 1987). For d i> 1 Gpc a cor- 



rection for cosmological curvature is necessary. The exact 
waveforms for eccentric orbits and resultant detectability of 
the radiation are more complicated (see for example Cutler 
and Flannagan 1995, Ryan 1995, Shibata 1994, Junker and 
Schafer 1992). 

What, then is the likely rate of such events? Hils and 
Bender (1995) recently considered the rate of capture in 
M32-like systems with isothermal cusps due to the diffu- 
sion of compact stars into the black hole, the diffusion being 
dominated interaction with main sequence stars on more 
loosely bound orbits. This work extends the consideration 
to non-isothermal systems and includes the contribution of 
both diffusion and large-angle scattering to the capture rate, 
and the mutual interaction of compact stars in tightly bound 
orbits. 

2 CUSPS AND SWALLOWING RATES 

Consider a massive black hole in the center of some galaxy. 
In general, the black hole will be surrounded by some 
(evolved) stellar population with a range of masses corre- 
sponding to some initial mass function and an associated 
population of evolved remnants stars, white dwarfs, neu- 
tron stars and low mass black holes. The stars will have 
some underlying density profile, with the density p(r) typi- 
cally well approximated by a (broken) power law (Kormendy 
and Richstone 1995) with some "break radius", rj,. For the 
galaxies which are the best candidates for harbouring black 
holes of mass ~ Mq, r(, is observed to be a few parsecs. 
We assume for now that the stellar distribution is spherical 
and the velocity distribution isotropic. 

We define a radius of influence of the black hole, = 
GM^/fTc, where a c is the one dimensional dispersion for 
f ^ r h ■ For most systems of interest the dispersion is ap- 
proximately constant for radii larger than but comparable 
to r/j , and is well defined. For the systems we will be 
interested in, ~ Mg and a c ~ 150 km s -1 . Defining 

&166 = fc/166kms -1 , it is useful to scale 



"166 

We define a characteristic dynamical time scale, t^y-a = 
r^/cTc, the time scale for a star with characteristic 1-D 
velocity to cross ; for systems considered here t^y-a = 
6000(r/j /pc)/<ri66) years. Inside the density profile is 
modified by the presence of the black hole (Peebles 1972, 
Bahcall and Wolf 1976, 1977, Cohn and Kulsrud 1978, 
Young 1980, Shapiro 1985, Murphy et al. 1991, Quinlan 
et al. 1995). 

A stellar population "relaxes" on some characteristic 
time scale, tf> = i/2<Tc /nG' 2 m*p\og(l/2N*), where N* is 
the number of stars (see eg. Binney and Tremaine 1987). 
Scaling to a normalised density pg = p/10 e pc -3 , we find 
tf> = 10 10 cr 3 gg//96 years. Note that tf> ^ tff(= 1-5 x 
10 10 years) is possible for reasonable values of <r c , p. 

If the lifetime of the stellar population in the center 
is longer than few times t/j, then the stellar population is 
relaxed; in particular mass segregation occurs due to dy- 
namical friction. The lifetime of the population of stars at 
the center of the galaxy is obviously ^ tjj, indeed it may be 
■C tjj if there has been recent star formation with associ- 
ated mass loss, or if the black hole arrived or formed in the 
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center of the galaxy recently. 
2.1 Cusps 

A self-consistent cusp of stars forms around a central black 
hole with some density profile p(r < r/j) oc r~ C (Peebles 
1972, Bahcall and Wolf 1976, 1977, Young 1980, Shapiro 
1985, Quinlan et al. 1995). In general C is a function of r, 
with different processes modifying the local slope at different 
radii and densities. The black hole may also induce a small 
anisotropy in the velocity distribution (Goodman and Bin- 
ney 1984, Quinlan et al. 1995). The cusp extends to some 
inner radius, r jn , where the star are effectively destroyed or 
swallowed. 

It is useful to parametrise the cusp slope as C(r) = 
3/2 +p(r). With t R oc a(rf p~ X inside r h , er(r) oc r -1 / 2 
and hence tR oc r p . A black hole growing adiabatically 
in a flat isothermal core, with tfj ^> tjj, induces a cusp 
with p(r < r h ) = (Young 1980, Quinlan 1995). This has 
the interesting property that t/j is independent of radius 
inside r^. For physical cores C < 2.5 and hence p < 1 
for plausible cusps whether relaxed or not (see Quinlan et 
al. 1995, Sigurdsson et al. 1995, for discussion). For a 
relaxed population of equal mass stars in the central region 
(tR( r h) "C min{tff , tj orm }), p(r) = 1/4 and the relaxation 
time decreases slowly with r, while for a relaxed cusp with 
a range of stellar masses, p ~ 0.3, with steeper cusps for the 
more massive stars (Bahcall and Wolf 1976, 1977, Murphy 
et al. 1991). 

Two further quantities are fundamental in determin- 
ing the structure of the cusp around the black hole. The 
collision radius, r co u = 7 x 10 16 M6 cm for main sequence 
stars (smaller by ryyD l r * f° r white dwarfs, effectively zero 
for our purposes for neutron stars and stellar mass black 
holes), is the radius at which star-star encounters cannot 
lead to large angle elastic scattering as the characteristic 
encounter velocities exceed the surface escape velocities of 
the stars (Frank and Rees 1976). Inside r co u two body re- 
laxation is by definition ineffective and t/j becomes large 
for the relevant population. Note that r co u ^> rj< for both 
stars and white dwarfs for the black hole masses of interest 
here. Another quantity is the black hole wandering radius, 
r w , which measures the root mean square displacement of 
the black hole from the galaxy's center of mass due to the 
discreteness of the potential of the galaxy. For a black hole 
in an isothermal stellar core r w ftj \J m* /M^r^ (Bahcall and 
Wolf 1976). For other core profiles the wandering radius 
is less well defined (Quinlan private communication). As- 
suming the graininess due to the stellar mass distribution 
dominates r w , then for the density profiles and black hole 
masses of interest here, ^> r w ^> ry. Consequently, de- 
pletion of the main sequence stellar population due to tidal 
disruption is not limited by loss-cone diffusion into rj< as the 
black hole wanders on dynamical time scales short compared 
to tfj. The short period random walk of the central black 
hole does not affect our rate estimates because the rate es- 
timates are stochastic, averaged over time and the ensemble 
of black holes. The short time scale fluctuations that move 
the black hole away from a compact remnant star that would 
otherwise have entered the gravitational wave loss-cone, are 
balanced on average by central black holes random walking 
into the path of a compact remnant star that would other- 



wise not have entered the gravitational radiation loss-cone. 

Another critical radius is set by the radius at which 
mergers of main sequence stars with each other become im- 
portant, r m . The rate for star-star mergers per star is sim- 
ply 

R cc = 7rr ' P(r) " (r) , (5) 

for r m < r/j, assuming a i> v esc , where v esc is the es- 
cape velocity at the surface of the star (true for r <C r^). 
For white dwarf-white dwarf and white dwarf-neutron star 
merger gravitational focusing increases R C c by a factor 
(1 + 2Gm* /r*<r(r) 2 ). In our units, 

fl^W (£)'(£)(?)"""' P.*. ,„ 

By definition, depletion of stars through merger becomes 
important when R C c ~ t]^ 1 , provided tf> i> tjj, else deple- 
tion only becomes effective at r co u . For solar mass main 
sequence stars, r m ~ 10 -1 r/j ftj r C oll i as expected. In- 
side min{r m , r co u } the cusp induced by the black hole is 
flattened by the depletion of main sequence stars due to 
mergers, with the main sequence density profile flattening 

to p(r < r m ) oc r~ D , where D ~ - 1/2 (p = -1 3/2). 

(Lightman and Shapiro 1977, Murphy et al. 1991, Rauch 
1995). 

While the main sequence stellar density inside r m is 
flat, this radius has no special significance for the compact 
remnants which would maintain a steep profile inside r m 
(see eg Murphy et al. 1991), with 

Pcr(r < r m ) = fcrp(r m ) > ( 7 ) 

where f cr is the fractional density of white dwarfs (or neu- 
tron stars or stellar mass black holes) at r m , corrected for 
mass segregation if applicable. Compact remnants stars also 
merge with main sequence stars, as well as each other, the 
rate for merger with main sequence stars is smaller than the 
star-star merger rate by a factor f cr , and simulations of such 
high speed mergers suggest in most cases the compact star 
will remain after the merger (Shara and Regev 1986, Benz et 
al. 1989, Ruffert 1992). Inside r coll (MS) the two body re- 
laxation time scale for white dwarfs increases sharply as the 
effective density of bodies available for large angle scattering 
decreases by f cr . 

2.2 Dynamics of compact stars in the cusp 

2.2.1 The influence of the central black hole 

Consider now the dynamics of the stars in the cusp. The 
stellar orbits are dominated by the gravity of the central 
black hole, but are perturbed by the mutual interactions of 
the stars. Stars which venture too close to the black hole 
are either tidally disrupted, in the case of main sequence 
stars, or swallowed whole, in the case of compact stellar 
remnants. In the latter case, the flux of compact remnants 
into the black hole is due to stars in the cusp scattered into 
orbits that either plunge them directly into the black hole, or 
such that the peribothron is small enough for gravitational 
radiation to shrink the orbit more rapidly than interactions 
with other objects either scatter the star away from the black 
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hole again, or put it on an orbit that plunges straight into 
the black hole. 

To first non-zero order, gravitational radiation leads to 
a decrease in orbital energy, E, and angular momentum, L, 
with 



dE_ 
dt 

and 

dL_ 

dt 
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where r p = a(l — e) is the peribothron of an orbit with semi- 



=2)3/2 

37/96e 4 )/(l + e) 5 and g' (e) = (1 + 7/8e 2 )(l ■ 



major axis a, eccentricity e, /'(e) = (1- 



(l + 73/24e 2 + 
■e) 3 /(l + e) 2 



Note that /' is less sensitive than g 1 to the eccentricity for 
e ~ 1 while both E and L are equally sensitive to r p . 

Frank and Rees (1976) considered the diffusion of main 
sequence stars in cusps around central black holes in galaxies 
with isothermal cores. We follow their argument for the scat- 
tering of compact stellar remnants for the full range of adi- 
abatic and relaxed stellar cusps around central black holes. 
We consider stars inside the cusp at radii r$ <C r r^. In 
order for the stars to be swallowed by the central black hole, 
their orbits must be within some critical radius r c (r). It is 
useful to define a "loss-cone", 8(r) = ^/2r mj - n /3r, where 
r min is the peribothron distance for the star at r. As the 
stars orbit about the black hole, the orbits are scattered by 
the inhomogenous potential they move in; it is useful to con- 
sider two regimes; where the scattering angle is small com- 
pared to 8, which we refer to as "diffusion", and scattering 
where the scattering angle is large compared to 8, which we 
refer to as "kicks" . The scattering in the respective regimes 
can be thought of as being due to the "Poisson noise" in the 
potential due to the discrete number of stars for diffusion, 
and as elastic scattering off individual stars in the case of 
"kicks" . In the absence of a central black hole, the scattering 
in and out of the loss-cone would be symmetric, with the 
flux into the loss-cone balanced by the flux out. In the pres- 
ence of a black hole there is an additional source of orbital 
evolution: the secular decay of the low angular momentum 
orbits due to gravitational radiation. Hence there is a net 
loss of stars to the black hole. 

The ratios of the "scattering time scale", t sca t and the 
time scale for decay through gravitational radiation, tow ■> 
t or (,(r), can be written 



to the orbital time scale 
2 



Mh 



where N*(r) oc f 



3/2-p(r) 



(10) 



is the number of stars interior to 
r. The time scale for decay through gravitational radiation 
is set by the energy radiation rate, given by 



taw 
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where a is the semi-major axis of the stars orbit about the 
black hole, and /(e) = (l + 73/24e 2 +37/96e 4 )(l+e)" 7/2 (l- 
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For the highly radial orbits we're interested in, 
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tscat is the time scale for scattering out due to "kicks"; 
the corresponding diffusion time scale is shorter by a fac- 
tor r m { n /r = 3# 2 (r)/2. In order for a star to be scattered 
to a small enough peribothron that it will be swallowed due 
to gravitational radiation, we require t sca t > tQfy, or that 
8 < 9 cr it- Hence for "kicks", we find 
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and for diffusion 
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For radii such that 9 cr n(r) > 1 the power-law cusp of stars 
is no longer present as stars are swallowed by the central 
black hole. 

The rate, R s , at which the stars are swallowed, either by 
scattering straight into the black hole, or by gradual shrink- 
age of their orbits is given by 



R s 



tscat 



(15) 



and we can solve for R s for a choice of cusp parameters and 
sum over the galaxy and black hole mass function for an 
estimate of the total rate. Substituting for 8 cr n 



2-2/13 , x -5//3 
m* \ r \ 1 



rs 



(16) 



P' 



where f3 = 5 for stars diffusing in the loss-cone, and f3 = 7 
for stars undergoing large angle scattering out of the loss- 
cone, and C 2//3 = 3/2(85tt/24V2) 2/7 « 2.7 for /3 = 7, and 
Cl ll} = 3/2(857r/24V2) 2 / 5 « 3.4 for (3 = 5. 



2.2.2 Single power law cusps 

If the cusp density is a constant power law, the number of 
stars, iV*(r), is given simply by 



N^r) = 4tt / p(r')r' 2 dr'; 



normalising to p(r/j) = Mq pc we find 



(17) 



1.2 x 10 7 p 6 f r h \ 3 f r \ 3/2 P 
Mr)= 3/2 -p m */M Q {—c) UJ • (18) 

Assume we have a single power law cusp, that p is inde- 
pendent of r for r < f5. It is useful to define r cr n, where 
S C rit{r C rit) = 1, then 



3-2p 

7—1 I — I r S- 



(19) 



3/2 -p 
1.2 x 10 7 p 6 J \r s 



A plot of r cr n vs. p is shown in figure 1, for M§, pe = 1, and 
r/j = 0.5, 1.0 and 5.0 pc respectively. The curves are shown 
for f3 = 7; r cr n is 25% larger for f3 = 5. 
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Table 1. Scaling of scattering rates. Assuming p(r^) = 
lO 6 M pc _3 , M h = M 6 and a = 166 km s -1 , r coll /r s = 2000, 
Rs(r) cx (r/rsP, 7 = (3 - 4p)/2 - (8 - 2p)/f3. 



Figure 1. r cr ^ vs p for single power law cusps for different . 

Then r crit « 600r s for p = 0, and r c „'t « r s for p = 1, 
much less than r m . Large angle scattering becomes ineffec- 
tive for white dwarfs at r co u PS 2000rs > r cr n, but neutron 
stars and stellar mass black holes may undergo large angle 
scattering to r ~ r$ • 

We consider the rate of capture for f3 = 5, 7 separately, 
still assuming a single power law cusp composed purely of 
compact remnants. Substituting in equation 17, we can 
write 
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Table 1 shows the resultant r cr n and capture rates for 
different slope cusps. 

The diffusion rate decreases with r for all p. The large 
angle scattering rate increases with r for p < 5/24. How- 
ever r cr n is much smaller than r co u for white dwarfs for the 
steeper cusps while for flat cusps the diffusion rate is rela- 
tively insensitive to r/r cr n. Note the rates shown in the last 
column are not the true rates as r cr n < r co u in all cases and 
large angle scattering is ineffective at these radii for white 
dwarfs. Figure 2 shows the capture rates for single power 
law cusps composed purely of compact remnants, for some 
different p, neglecting r co u . 

The fraction of stars that plunges directly into the black 
hole without gradual inspiral through gravitational radia- 
tion is R s (j3 = 5)/(R s (j3 = 5) + R s (j3 = 7)). For single 
power law cusps, this can be a small fraction; for more re- 
alistic cusps we expect 1/3-1/2 the stars to plunge rapidly 



Figure 2. Plot of R s (r) for different p, assuming single power 
law cusps composed entirely of compact remnants and ignoring 
existence of r co n . The solid line is for (3 = 7 (large angle scatter- 
ing) and the dotted line is for (3 = 5 (diffusion). 

into the black hole, with the remainder undergoing a more 
gradual inspiral with the orbit eccentricity decreasing. Most 
stars plunge from r ~ r m ^> r cr n and enter r < r cr n with 
eccentricities of ~ 0.999. 

For detection of gravitational radiation, the time the 
compact objects spend at orbital periods P ^ Pa{= P/10 4 s) 
is critical. The time to decay for highly eccentric orbits, 
using the quadrupole approximation, is well approximated 
by 



T GW = 10° (1 
for M 6 



„8/3 

Pa y 



(24) 



1, m* = Mq, note also tq\\? m. mj 
(see eg. Rajagopal and Romani, 1995). For e = 



Ml 



2/3 



-1 



, TO* 

0.99 at 

r = 2000rs the time to inspiral is ~ lO 4 years, for e = 0.999 
and r = 4000rs the time to inspiral is only 30 years and the 
orbit has no time to circularise substantially before reaching 
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periods P4. For steep (large p) cusps, where the rate is 
dominated by stars near r co u , diffusion will random walk 
e to somewhat lower values for a significant fraction of the 
stars, as discussed above. The eccentricity, ej, as the star 
approaches capture is smaller with, (f — ej) ~ 4 x (1 — e;), 
where e; is the initial orbital eccentricity after scattering, 
and the system may spend times of order fO 3 years at pe- 
riods ~ P4. As discussed in a recent paper by Rauch and 
Tremaine (f996), resonant effects may enhance the relax- 
ation rate of angular momentum in Keplerian potentials; 
however, this does not appear to be effective for the orbital 
parameters of concern here, because the relativistic preces- 
sion of the relevant orbits (with small peribothrons) is fast 
enough to destroy the resonances. Eccentricity evolution 
inside the loss-cone is therefore small. 

Stellar mass black holes may contribute strongly to the 
total observable rate. They spiral in an order of magnitude 
faster than white dwarfs and neutron stars and are observ- 
able for correspondingly shorter times at these frequencies. 
On the other hand, h c is an order of magnitude larger for the 
low mass black holes than for white dwarfs, and the volume 
for detection is ~ 10 3 times larger. If stellar mass black holes 
form in population II in significant numbers, (Jbh ^ 10 -4 ), 
they will provide a strong characteristic signal for LISA (see 
also Polnarev & Rees 1994). 

In practice cusps do not have simple single power law 
density profiles as we discussed above, nor are they com- 
posed solely from compact remnants. We now consider the 
expected capture rates in real cusps. 

3 REAL GALAXIES 

There are two classes of galaxies whose centres are good 
candidates for harbouring central black holes in the appro- 
priate mass range, which may be capturing compact ob- 
jects at an interesting rate: nucleated spiral bulges, such 
as that of our own galaxy; and the cores of compact dwarf 
ellipticals like M32. Both have low dispersion, steep high 
density central cusps, and probably contain central black 
holes of ~ 10 6 M . M32 was discussed in detail by Hils and 
Bender (1995), assuming diffusion dominated the capture 
rate and that the central cusp was unrelaxed and isother- 
mal (p = 0) It is likely that the core of M32 is in fact re- 
laxed, with p = 0.25, although observation at radii ~ 0, lr/j 
are consistent with p = 0. Systems like M32 are particu- 
larly promising sources because of their large central densi- 
ties and low central dispersion. We infer a capture rate of 
fa 3 x 10 -8 per year, dominated by large angle scattering 
from r fa 5 x 10 3 rs. A somewhat smaller diffusion/collision 
rate of 1.8 x 10 -8 was estimated by Hils and Bender (1995). 
Unfortunately the space density of dwarf galaxies like M32 
appears to be very small, probably as low as 10 -5 Mpc -3 
(Kormendy private communication, see also Gebhardt et al. 
1996), the total capture rate out to 1 Gpc is thus only about 

10-V 1 . 

3.1 Bulges of spirals 

The space density of nucleated spirals with the appropriate 
bulge mass is higher ~ 10 -2 ' 5 Mpc -3 , but the dynamics 
of their central regions are more complicated. At r/j typ- 
ical density profiles correspond to p Si 0.5 — 0.8. For the 



densities and dispersions seen in nuclei of spirals, the cusp 
of main sequence stars is flattened due to stellar mergers 
at r m Si 0.1r/j as discussed above. Inside r m the main 
sequence density profile is flat, PMs( r ) = / 9 o( r / r m) -1/ ' 2 , 
Po = P6( r h/ r m) 3 ^^~ p Si lO 8 M0pc -3 . The relaxation time 
at r m , is tn(r m ) ~ 3 x 10 9 years. The total number of main 
sequence stars inside r m is thus Nms sj 5 x 10 5 . The white 
dwarf-main sequence merger rate is smaller by a factor f C r , 
and it is likely at the encounter velocities seen in these cores 
that a white dwarf would emerge relatively unscathed from 
any such merger. Thus we expect the white dwarf (and neu- 
tron star and low mass black hole) density profile to remain 
steep. 

Including gravitational focusing, the cross-section for 
WD-WD mergers is R CC (WD) = 6 x 10" 5 (r h /r)- 2 - p t' 1 . 
Solving for r m (W D), and requiring R C c(WD) = 5 as white 
dwarfs interior to r m (MS) are replenished on tf{(r m ), we 
find r m (WD) « I0~ 2 r h < r co u(MS). The relaxation 
time scale for white dwarfs increases around r co u (MS) as 
WD-MS scattering becomes ineffective for two body relax- 
ation, but then decreases like r p inside r co n(MS), with 
r n(WD) <J tjj at r m (W D). Consequently, relaxation main- 
tains the white dwarf density profile at p ~ 0.3 down to 
r C oll(WD) ~ 5 x W- 3 r h ~ 2 x 10 3 r s . Inside r coll (WD) we 
expect the density profile of the white dwarf population to 
flatten out to pwD(r) oc (r/r m (WD))' 11 " 2 . 

The number of white dwarfs interior to r m (MS) is ap- 
proximately N\y£) = 10 4 . The density of neutron stars and 
low mass black hole does not flatten due to mergers, but will 
level off due to relaxation inside r m (WD) to p = 1/4 profile. 
Note that with the density profile flattened the relaxation 
time no longer decreases inside r co u (WD) Inside r co u (WD) 
we thus find t/j ~ tjj again and a strong (p i> 0.5) neutron 
star and low mass black hole cusp may persist all the way 
to Tcrit- 

The density profile of the neutron star and low mass 
black hole population follows that of the main sequence 
stars to r m (MS), p NS /BH = /jV S / BH Pe(r Aft ) ( ~ 3/2 ~ p) . At 
r m (MS) the profile may flatten to p = 0.25 — 0.3 as relax- 
ation becomes effective. The larger value is appropriate if 
the main sequence population is evolved; the lower value is 
appropriate for a younger main sequence population (as ap- 
pears to be the case in our own galaxy). The relaxed cusp 
profile persists to r co u (WD) at which point the tfj is large 
and the profile may steepen to p ~ 0.5 again. 

We calculated numerical large angle scattering rates and 
diffusion rates for piecewise power law cusps, where we al- 
lowed for the changes in the density profile inside r/j due to 
stellar mergers, changes in local relaxation time scales, and 
ineffectiveness of collisional relaxation inside r co u . Figure 
3 shows the capture rates inferred as a function of radius 
for some characteristic profiles expected in the centres of 
nucleated spirals. Typical inferred WD-MBH capture rate, 
through large angle scattering from r ~ 2r co u(WD) is 10 -7 
per year, for nucleated cusps with structural parameters like 
the Milky Way. Large angle scattering dominates the cap- 
ture rate, and the total rate is dominated by white dwarfs 
at radii ~ r m (WD) ~ 3 x 10 3 r s . With ~ 10 4 WDs inside 
r m the white dwarf population is not significantly depleted 
by mergers on time scales of tjj, and can be replenished by 
relaxation or stellar evolution as implicitly assumed in the 



Capture of compact stars by black holes 7 



Figure 3. The rate of capture of white dwarfs by a 10 Mq 
central black hole in a canonical nucleated spiral galaxy, as a 
function of radius. The rate peaks strongly at r ~ 2r co n(WD) 
at ~ 10 — 7 per year. 

derivation. 

If the space density of these bulges is 0.003 Mpc -3 , we 
expect ~ 1 captures per year within 1 Gpc. At any one time 
we may expect 10 2 systems in the process of capture within 1 
Gpc, or 0.1 — 1 system within 100 Mpc. Such systems would 
be easily detectable by LISA at 100 Mpc and might be de- 
tected near 10 -3 Hz out to 1 Gpc. The optimal signal is ex- 
pected from Mq PS 3: for smaller the gravitational wave 
amplitude becomes too small, for larger the frequencies 
become too low even for orbits close to r$ . Low central dis- 
persion leads to higher capture rates, favouring systems with 
intrinsic density cusps similar to Hernquist profiles, low cen- 
tral dispersion and p = 0.8 inside r/j. Whether such systems 
are prevalent in nature is an open question. 

If the IMF is flat and the neutron star fraction is higher 
than the 0.2% expected from a Salpeter slope mass func- 
tion, then neutron star captures may be competitive with 
the white dwarf capture rate due to the steep density pro- 
file expected for the neutron stars at all radii. If WD-WD 
mergers lead to accretion induced collapse and neutron star 
formation the expected capture rate may be dominated by 
neutron stars formed in the cusp. R S (NS) i> 10 -7 are possi- 
ble for /jvs ~ 1%i m which case the neutron stars contribute 
significantly to the total GW rate. However, neutron stars 
are expected to be born with natal kick velocities greater 
than (Ti66 an d the central potential may not be deep enough 
for many neutron stars to remain in the center after for- 
mation, depressing the effective f cr . Low mass black holes 
are detectable to larger distances, but presumably occur in 
much smaller numbers. LMBHs are, however, presumed not 
to receive natal kicks. 

If there is a substantial population of stellar mass black 
holes, the capture rate in spirals nuclei may be as high as 
10 -5 per year per galaxy for 10 Mq black holes, with a global 
number fraction of 2 x 10 -4 for the black holes (see figure 
4). With less than 10 3 black holes in the central regions, 
such a rate is clearly not sustainable, but would be possible 
following a nuclear star burst, for a period of 10 8 years or 
so. There is some evidence that the Milky Way underwent 
a nuclear starburst in the last 10 9 years, if this is typical 
of nucleated spirals, the population averaged rate for stellar 



mass black hole mergers within 1 Gpc might be as high as 
10 -6 per year per galaxy, in which case we might detect ~ 10 
such system coalescing per year with LISA. The expected 
initial eccentricities of the low mass black holes are ~ 0.999, 
and inspiral is rapid. 

3.1.1 Binaries, anisotropy and triaxiality 

There are likely to be some stellar binaries in galactic bulges. 
However, inside the density and dispersion are large and 
only binaries with semi-major axis ~ r* are hard enough 
not to be broken up by encounters with other stars. Such 
binaries cannot lead to a larger capture rate as the semi- 
major axis is <J r cr n and encounters with the binaries are 
no more effective in scattering them into the loss-cone than 
comparable single star scatterings. A significant fraction, /(,, 
of white dwarf-white dwarf binaries with semi-major axis 
~ Rq would lead to an enhanced WD-WD merger rate, by 
a ratio of a\YD-WD I RwD ~ IOO/5. 

The black hole also polarizes the stellar distribution 
inducing a finite tangential anisotropy at small radii (see 
Quinlan et al. 1995 for discussion). There is consequently 
some bias towards circular orbits in the vicinity of the 
black hole which will lower the estimated capture rate by 
~ 10%. As other processes can induce comparably mild ra- 
dial anisotropies (in particular ejection of single stars from 
the inner cusp by low mass black holes bound transiently to 
the central black hole), this does not change our estimate of 
the capture rate. 

In the inner cusp the stellar distribution is forced to- 
wards sphericity by the black hole potential (there might in 
general be some modest rotational flattening). On larger 
scales the spheroid is most likely triaxial, possibly strongly 
triaxial. Orbit diffusion in triaxial potentials may be a 
strong factor in replenishing the loss-cones of central mas- 
sive black holes; this effect is irrelevant, as the black hole 
mass is low enough that its "wandering" ensures the loss- 
cone remains filled independent of the shape of the stellar 
distribution. 

As discussed by, for instance, Syer, Clarke and Rees 
(1991) a main sequence star on an eccentric orbit with peri- 
bothron somewhat larger than rj< can cumulatively lose en- 
ergy by impact on an accretion disk, so that it ends up on 
a tightly bound circular orbit. This process depends on the 
geometric cross section of the star and would be correspond- 
ingly less effective for a compact object. The presence of a 
modest mass gas accretion disk around the central black hole 
will therefor not affect our present estimates of the capture 
rate for compact objects. 

3.2 Stellar mergers 

The rate of WD-WD mergers is high within the black hole 
induced cusp. The merger rate per white dwarf at r m (WD), 
Rcc(WD) = 3 x 10" 10 per year, by definition. There are 10 4 
white dwarfs interior to r m (WD) given the assumed cusp 
parameters. The integrated merger rate inside r m (WD), is 
N C (WD) ~ 10 -6 y -1 , enough to destroy every white dwarf 
in the cusp in a Hubble time. As the cusp is relaxed at 
r m (WD) this is not a concern and the white dwarf popula- 
tion can be replenished both through local stellar evolution 
and replenishment from outside r m (WD). The expected 
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Figure 4. The rate of capture of LMBHs by a 10 6 M central 
black hole in a canonical nucleated spiral galaxy, as a function 
of radius. The rate peaks strongly at r ~ lO^rg at ~ 10 — ^ per 
year. 

number of WD-WD mergers within 100 Mpc is 10 per year, 
with a NS-WD merger expected once every few years within 
the same volume, depending on /jvs- As noted above, a 
substantial fraction of WD-WD binaries could increase the 
WD-WD merger rate by an order of magnitude. 

If WD-WD mergers lead to accretion induced collapse, 
rather than total disruption of the white dwarfs, then WD- 
WD mergers may lead to a substantial neutron star popula- 
tion in the centres of galaxies, producing up to 10 4 neutron 
stars inside r m in a Hubble time. This would on on average 
double the estimated coalescence rates of white dwarfs and 
neutron stars with the central black hole. 

The white dwarf-main sequence merger rate is ~ 10 -5 
per year, and we assume a compact remnant remains after 
the collision. If the white dwarf is temporarily bound to the 
main sequence star after merger, a common envelope phase 
may ensue, producing a large, luminous stellar object with 
a white dwarf core. If such objects last 10 7 years, we would 
expect ~ 10 2 to be observed at any one time in the inner 
0.1 pc of a typical nucleated bulge. 



4 CONCLUSIONS 

We estimate the capture rate of compact stellar remnants 
by massive black holes in the centers of power law bulges of 
nucleated spirals and compact dwarf ellipticals. 

The capture rate in M32-like ellipticals is high, and 
likely to be dominated by large angle scattering by other 
stars in the tightly-bound cusp. The net rate for M32 like 
systems is likely to be of order 10 -8 per year, but the total 
rate in the galactic neighbourhood is probably small due to 
the dearth of such systems. 

The total rate is likely dominated by the power law 
bulges of "ordinary" spiral galaxies such as our own Milky 
Way. The total expected rate is somewhat lower than naive 
estimates due to the flattening of the white dwarf density 
profile through WD-WD mergers, with likely event rates of 
order 10 -8 per year per galaxy. However, due to their rel- 
atively high space density such systems dominate the total 
observable rate out to 1 Gpc. 



We conservatively estimate a minimum of 0.1 captures 
per year out to 1 Gpc, with perhaps 10 2 systems observable 
at low frequencies in the early stages of capture at any one 
time. Both the burst and periodic signals should be de- 
tectable by proposed gravitational radiation observatories 
such as LISA, out to few hundred Mpc, with a character- 
istic signal from the high eccentricity orbits. If low mass 
black holes are present in significant numbers, in the cen- 
ters of galaxies, then the signal from LMBHs captured onto 
the central black hole through large angle scattering may be 
an order of magnitude larger still, averaged over the local 
population of galaxies. 

The WD-MS merger rate is estimated at 10 -7 per year, 
with O(10 2 ) luminous stellar objects descended from such 
mergers observable in the inner nucleus at any one time. 
The WD-WD merger rate is estimated at 10 -6 per year per 
galaxy, with several per year expected in the local super- 
cluster. Such events may be detectable through X-ray or 
UV flaring by current space based observatories. Chemical 
contamination of the bulge from such events may also be 
significant (Khokhlov and Novikov in preparation). 
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Figure 1. r cr ^ vs p for single power law cusps for 
different . 

Figure 2. Plot of R s (r) for different p, assuming single 
power law cusps composed entirely of compact remnants and 
ignoring existence of r co n . The solid line is for (3 = 7 (large 
angle scattering) and the dotted line is for (3 = 5 (diffusion). 

Figure 3. The rate of capture of white dwarfs by a 
10 6 M central black hole in a canonical nucleated spiral 
galaxy, as a function of radius. The rate peaks strongly at 
r ~ 2r co n(WD) at ~ 10 — ^ per year. 

Figure 4. The rate of capture of LMBHs by a 10 6 M 
central black hole in a canonical nucleated spiral galaxy, as a 
function of radius. The rate peaks strongly at r ~ lO^rg at 
~ 10 — ^ per year. 
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